Charge dynamics in the Mott insulating phase of the ionic Hubbard model 
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I extend to charge and bond operators the transformation that maps the ionic Hubbard model at half 
filling onto an effective spin Hamiltonian. Using the transformed operators I calculate the amplitude 
of the charge density wave in different dimensions D. In fD, the charge-charge correlations at large 
distance d decay as d~^ In"'''''^ d, in spite of the finite charge gap, due to remaining charge-spin 
coupling. Bond-bond correlations decay as (— f )'*d~^ In"^''^ d as in the usual Hubbard model. 



I. INTRODUCTION 

The ionic Hubbard model (IHM) has been proposed 
in the 80 's for the description of the neutral-ionic tran- 
sition in mixed-stack charge-transfer organic crystals. 
In the 90 's the interest on the model increased due to its 
potential application to ferroelectric perovskites."^"^ The 
model in any bipartite lattice can be written as: 

H = Ho + Ht; Ho = jJ2^-iyn,, + uJ2n^]'rl,^ 

i<7 i 
Ht = -tY,c\^Sa^.,„^ (1) 

No- 
where i-\-5 denote the nearest neighbors of site i, and odd 
i correspond to the sublattice with on-site energy — A/2. 

At zero temperature, in the strong coupling limit i = 0, 
the system is in the Mott insulating (MI) phase for 
[/ > A (all sites are singly occupied), and in the bond 
insulating (BI) phase for J7 < A (all sites with odd i 
doubly occupied). In one dimension (ID), for finite i, 
field theory calculations pointed out the existence of a 
third intermediate phase, the spontaneous dimerized in- 
sulator (SDI) between the other two.* For fixed A, as 
U increases, first a transition at Uc (involving mainly 
charge degrees of freedom) between the BI and the SDI 
takes place, and then for Us > Uc, the spin gap closes as 
in usual Kosterlitz-Thouless transitions. The physics has 
some similarities to that of the extended Hubbard model 
for small U and nearest-neighbor repulsion V.^~^^ 

The phase diagram of the ID IHM has been calculated 
numerically using the method of jumps of Berry phases,^ 
which in this case coincides with the method of cross- 
ings of appropriately chosen excited levels.®'"'^'' The charge 
Berry phase is a measure of the polarization and there- 
fore detects the ferroelectric charge transition, "'^^ while 
the spin Berry phase jumps at the point where the spin 
gap closes. For small t, Ug— Uc = 0-Qt was obtained for 
the width of the SDI phase. ^ 

Due to the fact that the spin gap is exponentially small 
for U smaller but near Us, direct numerical calculation 
of it or of appropriate correlation functions are unable 
to find a sharp transition at Us-^^''^^ The same difficulty 
happens in the Hubbard chain with correlated hopping, 



where the existence and position of the transition is con- 
firmed by field theory^^'^^ and exact^* results. In the 
IHM some controversy remains about the existence of Us 
and the nature of the phase for U > Uc in the case of 
absence of a second transition. Wilkens and Martin-'^^ 
suggested that the MI phase does not exist (the spin gap 
remains open for any finite U). This is in contradiction 
with the strong coupling expansion for < <C ?7 — A, which 
maps the IHM onto an effective spin Hamiltonian H with 
closed spin gap.^-^"* Wilkens and Martin^'^ argued that 
since H has an additional symmetry with respect to H 
(translation to a nearest neighbor Tg), H seems to have 
lost part of the physics of H , which might be essential. 

Recently, a detailed study of different observables of 
the ID IHM was made, using density-matrix renor- 
malizarion group (DMRG) , including a careful finite-size 
scaling analysis. The results indicate the present of two 
transitions and are consistent with the phase diagram 
found using topological transitions,^ although a smaller 
width of the SDI phase is suggested. Surprisingly, the 
authors find a power law decay of charge-charge correla- 
tion functions in the MI phase, in spite of a the presence 
of a charge gap ~ J7 — A. 

Here I report calculations of charge expectations val- 
ues, and in one dimension, the long distance behavior of 
charge-charge and bond-bond correlation functions using 
the mapping to a spin Hamiltonian, valid for t ^ ?7 — A. 
This might seem surprising at first sight, since charge 
fluctuations are frozen in H . However, the relevant infor- 
mation is contained in the transformed operators. This 
has some similarities with the case of the cuprates, in 
which Cu and O contributions to the photoemission spec- 
tra were successfully calculated using effective one-band 
models. Also, while at — J model is enough to ex- 
plain the observed dispersion in photoemission measure- 
ments in Sr2Cu02Cl2, the mapping of the operators^^'^^ 
is crucial to explain the observed intensities.^^ 

The transformation on the operators is performed in 
Section II. Section HI contains the results for expecta- 
tion values of different quantities, and Section IV is a 
discussion. 



1 



II. THE CANONICAL TRANSFORMATION 

In this Section, I use the canonical transformation plus 

projection onto the low-energy subspace that maps H 
into a spin Hamiltonian H, to transform charge and bond 
operators in one direction (i5 = ±1) 



rii = + Tiii, 

=5I'^(^^+<5<TC^<T + H.C.), 
6cr 



(2) 



and discuss some symmetry properties. For our purposes, 
it is enough to work up to second order in t/{U — A). 
Thus 

H = Pe-^He^P = P{H + [H, S] + ^[[H, S],S] + ...)P, 

(3) 

where P is the projector over the ground state subspace 
of Ho, in which rij = 1 for all i, and S is chosen in 
such a way that linear terms in the hopping term Ht are 
eliminated: Hi + [H^, S] = 0. Taking matrix elements of 
this equation between eigenstates of Hq, one obtains: 



{n\S\m) 



{n\Ht\m) 

Efa En 



(4) 



where Ej is the energy of the cigcnstate |j) of Hq. From 
Eqs. (3) and (4), proceeding in a similar way as done 
below, one obtains the known result:^ 



H = 



([72 _ A2) 



(5) 



{ij) 



The most important correction to H in higher order is a 

next-ncarcst-ncighbor exchange of order t^, which does 
not affect the physics for f <C [/ — A.^'^^ 

The effective Hamiltonian H to all orders in t is in- 
variant under a nearest-neighbor translation Tg, while H 
is not. 2 This is simply a consequence of the fact that 
to all orders H \s a purely spin Hamiltonian, since the 
charge degrees of freedom are frozen (n^ = 1 for all i). 
Then H is invariant under the electron-hole transforma- 
tion Teh ■ c^iu cFCia, which leaves invariant all spin 
operators. Since the original Hamiltonian H is invariant 
under the product TsTeh, then H should also be invari- 
ant under Ts- However other transformed operators, like 
fii below, are not invariant under Tg. 

The transformed charge operator is: 



Pe-^n,e^P ^ Pin, + [ni,S] + ^[[m, S], S])P 



1 + P{SS - SniS)P. 



(6) 



In the second equality I used the fact that S applied 
to any state of the ground state manifold of Hq, gives an 
excited state. The only terms of —PSriiSP which do not 



cancel with PSSP are those in which = or rij = 2 
after the first application of S. Using Eq. (4): 



5(T(t' 



t 



4 



u + i- 



-1)'a) 4-|-5o-^»o'4<7'^»+5o-']-P- 



(7) 



The first (negative) term correspond to charge rii = 
in the intermediate state and the excitation energy is 

Er,-Em = U+ (-1)*+^ A/2 - (-l)*A/2 = U- {-if A. 
The sums over spins can be transformed as follows: 

= J2 PnjAl - n,^)P - S+S- - 575+ = 



— — 2Sj • S,'. 



And replacing this in Eq. (7) : 
2[/A<2 



rii 



1-i-iy 



(C/2 - A2)2 



^(1 - 4Si ■ Si+s). (9) 



This equation relates a charge operator with a purely spin 
operator. Although this might seem surprising at first 
sight, it is clear that in the strong coupling limit, virtual 
charge fiuctuations are inhibited only for a ferromagnetic 
alignment of the spins. This shows that for finite U, some 
charge degrees of freedom remain coupled with spin, in 
contrast to the usual Hubbard model (A = 0). 

Using Eq. (2), to leading order in t, the transformed 
bond operator is: 

bi=Pe-^bie^P^P[bi,S]P 

+ U _ ^_-^y^ 4a'(^i+S'^'4+5a<^i'^ + H.C.]P. 

The first term correspond to acting first with S creating 
an excitation with a doubly occupied site at i, and an 
empty site dXi + 5, and then with 6j, restoring occupation 
1 at each site. Using Eq. (8), I obtain: 



-8tU 



C/2 - A2 



(10) 



The second member is a local measure of the asymmetry 
between the spin "bonds" involving site i in the direction 
of 6. Here there is no essential difference with the result 
for the ordinary Hubbard model. 



(8) 
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III. OBSERVABLES AND CORRELATION 
FUNCTIONS 

For i <C ?7 — A, using Eq. (9) the amplitude of the 
charge density wave is given by: 



= 2a^(l-4(Si-Si+5)^), 



(11) 



where the subscript in the expectation values indicates 
the Hamiltonian with which they are calculated and 



2UAt^ 

([/2_ A2)2- 



(12) 



In ID, from the Bcthc ansatz solution 1/4 — (S, • 
Si+s)H = ln2 ^ 0.69.26. This gives A ^ 11.09a. This re- 
sult has also been obtained using HcUman-Fcynman theo- 
rem and is in very good agreement with DMRG results. 
For the square lattice (Sj • Si+s) h — —0.3347,2'^ and then 
A = 18. 7o. In three dimensions, spin waves is a good 
aproximation and gives (S, • Sj+a)^ = —0.30,^^ what 
leads to A = 26.3a. 

The charge-charge correlation function is given by: 



Cd = {n^ni+d)H - {rii) H (ni+d) H ■ 



(13) 



To leading order, I can approximate {'nini+d)H — 
{nihi+d)H- The neglected term (see Eq. (6) ) 



(Pe-'^n,e'^(l - P)e- 



-dC P) involves charge fluc- 



tuations across the energy gap and therefore should 
lead to an exponentially decaying contribution, as in 
the usual Hubbard model. Using Eq. (9) in the form 
hr^l + (-l)HaSi • Si+8 - z{-iya + Oit"^), where z is 
the number of nearest neighbors, and dropping the sub- 
script H, one obtains to leading order in t: 



Cd = 16a\-lfJ2[iiSi ■ Si+5)(Si+d • Si+d+s)) 

SS' 

-(S,; ■ Sj+5)(Sj+d • Si+(j+5)]- (14) 

Note that the terms O(f^) in rij cancel exactly in Eq.(13). 
Thus, it is enogh to include terms up to order t"^ in fij to 
obtain the result up to order in Cd- Using symmetry, 
I can write: 



{{Si ■ Si+s){Si+d ■ Si+d+s)) — 3(5*0 S'fS'JS'l 



(15) 



and similarly for the other term inside square brackets of 
Eq. (14). 

In ID, the leading power-law decay can be determined 
using bosonization expressions for the XXZ model and 
correlation functions of the gaussian mo del. How- 
ever, at the isotropic point, there are important logarith- 
mic corrections. They can be calculated using expres- 
sions derived by Giamarchi and Schulz.^*^ Using S? w 



cos(7r^ -I- ^/2ip{x)) for the slowest decaying part of S'f , 
the operator product expansion SfSfj^-^ ~ cos(2-/2ip{x)) , 
and the results of section III B of Ref.^°, I obtain: 



{S^^SISlSl+s)^d-Hn-^d. 
Using symmetry again: 



(16) 



i^O^S^d^d+s) — gii^O^S + "^0 ^i)i^d^d+5 + ^d+s)) 

-{S^SISISI^s)- (17) 

The first term in this equation turns out to be the dom- 
inant one, and therefore I explain it in more detail. Per- 
forming a Jordan- Wigner transformation from spin op- 
erators to fcrmions with annihilation operators aj , going 
to the continuum limit using aj = P L[x) + {—iyR{x), 
with X = ja, and then bosonizing one gets: 

^{Spi+s + H.c.) ^ ^((a]a,+, + H.c.) 

6 6 

Y,\i^L^x)L{x + a6) + {-iYR^{x)R{x + aS) + 

5 

{-iyi-ifLl{x)Rix + aS) + {-l)H^ R\x)L{x + a6) + H.c] 

+<-"^«'fl+«- 
where c is a nonuniversal constant of the order of one. 
I would like to remark that the bosonic fields y, 6 are 

not pure spin operators of the original Hamiltonian Eq. 
(1), because they come from operators already dressed by 
the transformation. The expression of the transformed 
original fermions for A = and any filling to lowest order 
in t/U can be found in Refs.^^'^^. 

The first two terms of the last member of Eq. (18) lead 
also to contributions of order 1/d"^ (without including 
logarithmic corrections). Using^'^: 



(: cos(-\/2i^(x)) :: cos{V2ip{y)) :) 



1 



|a; — y\ h^^^ \x — y\ 
(19) 



together with Eqs. (14) to (19), I finally obtain (except 
for some factor of the order of one): 



7-/2 A 2j-4 

^^"'' (£/2-A2)4 ^''^"''^'^- 



(20) 



Note that factors (—I)'' cancel, and then in spite of the 
even-odd oscillations in (%■), these oscillations are absent 
in the charge-charge correlation functions. 

Concerning the bond-bond correlation functions, using 
Eqs. (10), (19) and E5 <5Si • ~ cos{V2ip{x)) , one 
obtains: 
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{bihi+d) « 64 



UH 



2+2 



(t/2 - A2)2 



{-ifd-^hi-^'^d. (21) 



These are much smaller than the spin-spin correlation 
functions, which behave like {SlS^j^^ w dr^lri^^^ d at 
large distances, and have a larger prefactor. 



IV. DISCUSSION 

I have calculated several quantities of the ionic Hub- 
bard model, using a mapping to a spin Hamiltonian valid 
for U > Is. and t <C J7 — A. In this limit the system is 
expected to be in the MI phase. This is predicted to be 
the case in ID,^'^ and in fact the absence of a spin gap 
in the effective transformed Hamiltonian H is consistent 
with this.2'1'* Also, the bond-bond correlation functions 
decay faster than the spin-spin ones at large distances, 
in contrast to what is expected in the SDL 

By construction, the symmetry of the observables is 
independent of the basis in which they are calculated, 
and this becomes clear after transformation of the oper- 
ators. For example, the charge operator is invariant 
under translation of one lattice parameter T^, while H 
is not. These symmetry properties are interchanged in 
the transformed operators fii and H , and the expecta- 
tion value {riiiH — is not invariant under Ts. As 
a consequence, for any finite U there is a charge den- 
sity wave. Its amplitude iov t <^ U — A is proportional 
to t'^AU/{U'^ - A2)2 but with a coefficient larger than 
20 which depends on dimension, and is near 53 in three 
dimensions. 

The charge-charge correlation functions in ID decay as 

d^^ In^'"*^^ d as a function of distance d. for large d. with 
a prefactor proportional to f^A^, in excelent agreement 
with recent DMRG calculations.^® This result for a finite 
charge gap, is in marked contrast to the general behavior 
found in ID systems and is a consequence of the fact that 
coupling between charge and spin dynamics remains at 
low energies in the original basis (effective dressed charge 
and spin fields are separated at low energies). This fact 
is at variancx; with the low-energy spin-charge separation 
that takes place in translationally invariant models with 
quite general nearest-neighbor interactions^'' (including 
the extended Hubbard model with nearest-neighbor re- 
pulsion and the Hubbard model with correlated 
hopping^"'^^'^''). An important difference is that in weak 
coupling the charge-spin interaction is a relevant opera- 
tor in the IHU}^ 
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